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ELECTROMAGNETIC  SCATTERING  PROPERTIES  OF  FINITE 
CYLINDERS  AND  SPHEROIDS 


L .  INTRODUCTION 

This  report  describes  the  work  done  and  the  progress  made  la  determining 
the  electromagnetic  (EM)  scattering  properties  of  3mnli  non-spherical 
particles.  In  particular,  the  purpose  of  the  contract  was  to  extend  and 
evaluate  the  modified  Shifrin  integro-dif ferenttal  aquation  solution!  to  the  EM 
scattering  by  finite  cylinders  and  spheroids.  These  solutions  are  then  used  in 
formulating  an  analytic  tool  for  the  investigation  of  multiple  scattering  by 
non-3phertcal  particles. 

The  solutions  for  scattering  by  finite  cylinders  and  spheroids  developed 
by  C.  Acqui3ta2  have  been  generalized  to  allow  for  arbitrary  polarization  of 
the  incident  EM  radiation,  polarization  scattered  directions^  and  arbitrary 
orientations  of  the  scattering  particles.  A  number  of  comparisons  with  "exact" 
theory^  have  been  developed  to  determine  the  range  of  applicability  of  the 
modified  Shifrin^  approach.  We  have,  in  addition,  developed  Internal 
electrostatic  solutions  foe  3hort  cylinders  so  that  the  Shifrin  method^  may  be 
used  to  compute  the  scattering  by  such  particles. 

A  very  detailed,  general,  double  scattering  problem  ha3  been  analyzed  and 
computations  performed  for  double  scattering  from  randomly  oriented.  Infinitely 
long  cylinders.  The  problem  of  double  scattering  from  finite  cylinders  and 
spheroids  has  been  analyzed;  however,  the  programming  of  the  solution  was  not 
completed. 

A  comparison  of  the  results  of  scattering  experiments  conducted  at 
Fairleigh  Dickinson  University  with  calculations  from  our  infinite  cylinder 
code  has  yielded  excellent  agreement^. 

2.  PERTURBATION  THEORY  FOR  SCATTERING  FROM  DIELECTRIC  SPHEROIDS  AND  SHORT 
CYLINDERS 

The  solution  of  Maxwell's  equation  for  the  scattering  of  a  plane  wave  from 
3n  Isotropic  and  non-magnet Lc  medium  can  be  expressed  by  the 
Integro-differential  equation 
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where  E(r)  i3  the  eLectria  field  far  from  the  scattering  medium,  c0  13  the 
propagation  vector  of  the  incident  wave,  and  m  is  the  Lndex  of  refraction  of 
tne  scattering  medium.  3y  expanding  the  VxVx  term  by  vector  identities  thl3 
equation  reduces  to 
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The  Shifrin-5  approach  to  sol/e  this  aquation  is  to  estimate  the  5  field 
inside  the  target  medium  and  use  thi3  as  the  start  of  an  iteractive  approach, 
tn  fact,  in  the  Rayleigh-Gans-Rocard  approximation  one  uses  the  external  plan 
wave  foe  the  internal  field.  Acquista's  extension^  of  Shiftin' 3  approximation 
begins  with  a  more  realistic  choice  for  the  lowest-orier  approximation  in  the 
internal  field,  namely,  the  internal  3oLution  for  an  infinite  dielectric 
cylinder  in  uniform  electrostatic  field.  We  adopt  thi3  general  approach  in  our 
work;  however,  we  first  need  to  modify  the  approach  for  the  case  when  the  light 
is  incident  on  the  cylinder  from  an  arbitrary  direction.  The  solution  inside 
the  cylinder  is  taken  to  be 
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where  Eeff  is  the  affective  polarizing  field  inside  the 
polarization  matrix  is 


cylinder. 


The  first-order  iterate  solution  is  found  to  be 
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for  the  electric  field  in  the  far  field.  The  derivation  of  thi3  equation  i3 
detailed  in  Reference  3.  Here  ko  is  the  propagation  vector  of  the  incident 
wave,  and  the  subscript  i.  means  only  the  portion  of  the  vector  perpendicular 
to  the  detector  direction  i3  to  be  taken.  The  pupiL  function  u  is  the  Fourier 
transform  of 
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where  a  is  the  radius  of  the  cylinder,  h  ia  the  length,  x*k  r-!c  ,  and  ||  anij. 
refer  to  the  cylinder  axis,  (see  Figure  1  of  Reference  3). 

The  second-order  approximation  is  found  to  be^ 
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where  the  /actor  F  is  defined  as 
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The  exact  determination  of  the  polarization  matrix  A  i3  found  by  solving 
the  electrostatic  problem  for  a  dielectric  finite  cylinder  in  a  constant 
electric  field.  This  problem  does  not  admit  a  closed  analytic  solution.  We 
have  3hown,  however,  that  by  using  the  infinite  cylinder  solution  as  our 
internal  trial  function,  the  iteration  procedure  converges  by  second  order  for 
aspect  ratios  ranging  from  infinity  to  approximately  23.  For  thi3  range,  the 
polarization  matrix  3eems  to  vary  only  slightly  as  is  shown  in  Reference  3. 


Rewriting  the  first-order  contribution  to  the  electric  field  we  get 
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A^  needs  to  be  expressed  in  the  detector  frame  (Figure  2  of  Reference  3)  and  we 
need  the  portion  of  A^  perpendicular  to  the  detector. 
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The  intensity  for  an  electric  field  polarized  in  both  the  initial 
and  final  state  is  defined  as 
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where  Egc^  2^  in  first  order,  and  l  is  a  unit  vector  in  the  direction  of  the 
final  3tate  of  polarization. 


Figure  3  of  Reference  3  shows  the  intensity  I, 


scattering 


infinite  cylinder.  The  3olid  curve  is  computed  from  first  order  and  the 
crosses(x)  are  the  result  of  the  exact  infinite  cylinder  computations.  It  is 
shown  in  Reference  3  that  this  excellent  agreement  is  a  consequence  of  the  fact 
that  Che  first-order  theory  becomes  analytically  equal  to  the  exact  Infinite 
cylinder  results  in  the  limit  of  the  height  of  the  cylinder  h  -*•  <»  . 

Since  our  goal  13  to  be  able  to  calculate  the  effect  of  multiple 
scattering,  we  investigated  the  range  of  cylinder  parameters  for  which  the 
first-order  solutions,  which  are  analytical  solutions,  converged.  Convergence 
was  assumed  when  the  second-order  iteration  differed  from  the  first  order  by 
lass  tnan  IX.  The  rationale  for  this  approach  i3  that  Che  first-order  solution 
requires  far  less  than  1  second  to  compute  an  entire  scattering  distribution, 
whereas  the  second-order  calculation  requires  about  20  minutes  from  the  same 
distribution  on  the  Per'xin  Elmer  3230.  In  multiple  scattering  where  the 
scattering  geometry  is  severely  complicated,  the  combination  of  the  geometry 
and  the  3acond-ordar  calculations  would  maxe  the  computations  impractical. 

The  results  of  our  investigation  show  that  the  first  order  is  within  IX  of 
tine  second  ordec  as  long  a3  the  aspect  ratio  is  20  or  greater,  and  the  phase 
shift  5  *  4m'im-l  ja/ A  is  I333  chan  2. 

The  scattering  cone,  characteristic  of  infinite  cylinders,  is  shown  to 
diffusa  as  tne  aspect  ratio  decreases  (see  Figures  4,  5,  5  and  7  of  Reference 
3). 


3.  SHORT  DIELECTRIC  CYLINDERS 

The  next  step  in  our  program  involved  the  study  of  short  cylinders,  with 
aspect  ratios  near  unity.  The  Shifrin  perturbation  theory3  relies  on  a 
reasonably  accurate  determination  of  the  electrostatic  polarization  matrix. 
Though  'xnown  for  the  infinite  cylinder  which  we  have  shown  wor’<3  reasonably 
wall  for  cylinders  down  to  aspect  ratios  of  20,  it  is  not  available  for  finite 
short  cylinders.  Na  have  U3ed  a  method  for  finding  4  by  a  technique  outlined 
in  a  boox  by  van  Sladel/ 


The  electrostatic  potential  <j>(rQ)  tn9ide  the  cylinder  i3  related  to  the 
surface  potential  as 
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where  <+>0  i3  the  ambient  potential,  the  integral  is  over  the  cylinder  surface, 
and  3/3n  is  the  outward  normal  gradient  to  the  surface.  It  follows  that  the 
electric  field  inside  the  cylinder  is 
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where  SQ  is  the  ambient  electrostatic  field.  This  internal  field  is  related  to 


the  electrostatic  polarization  matrix  as 
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and  we  are  led  to  the  following  equation  for  the  matrix  a  Laments: 
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where  a^E  =*  a^n  =»  a 2  and  a^  *  a^. 


Therefore,  the  matrix  A  is  known  if  we  know  the  surface  potential,  6  and 
it  is  determined  from  the  aquation 
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where  rg  is  on  the  cylindrical  surface.  We  solve  t'ni3  aquation  by  an  iteration 


process  3tartLng  with  $  a  33  the  first  guess.  Convergence  occurs  rapidly 
for  an  index  of  refraction  of  1.33  where  only  two  iterations  are  necessary. 


The  results  for  the  elements  a^,£  and  a^  are  given  in  Figure  7  of 
Reference  4  for  aspect  ratios  h/2a  *  1/2,  l  and  5.  The  values  are  given  along 
the  symmetry  axis,  but  they  are  similar  along  a  perpendicular  axi3. 
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It  is  observed  that  for  an  aspact  ratio  of  5,  the  elements  ara  quita 
constant  in  the  central  region  of  tha  cylinder  with  values  near  thosa  for  an 
infinite  cyiiadar.  The  values  of  the  elements  decrease  near  the  ands.  The 
same  general  variation  occurs  for  the  other  aspacts  ratios. 


4.  THE  atELECraiC  SPHENOID 

The  dialactric  spheroid  provides  U9  with  one  of  the  very  few  cases  for 
which  exact  solutions  axist?  and  thus  allows  us  to  avaluata  properly  the 
accuracy  of  the  Shifrin  approach^  and  the  ranga  of  convargence  of  the 
expansion.  For  the  spheroid,  the  pupil  function  a  is  well  known,  and  it  has 
tha  forms 

U(p)  =  V  f  ■[  a[p2  -  ]  |  prolate 

(15) 

U($)  *  V  f  |  a[p2  -  e2p(|  ]  |  oblate 

where  p^  and  p)(  ara  tha  components  of  p?  perpendicular  and  parallel  to  the 
symmetry  axis  of  the  targat,  a  is  tha  sami-major  axi3  of  ellipsoidal  cross 
section,  g  Is  the  eccentricity,  v  is  the  volume,  and 
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The  electrostatic  polarization  matrix  i3  also  well  known,  and  its  elements  for 
a  prolata  spheroid  ara 
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For  an  oblate  spheroid,  the  alamant3  ara 
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where  6  =  ^  .  Tha  3cattarlag  amplitudes  calculated  ara 
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The  rasalCs  ara  compared  with  the  exact  calculation  of  Reference  3  and 
appear  in  Reference  4  a3  Figures  2  through  4  for  the  respective  aspect  ratLos 
L/2,  1  (spnara),  and  5  taking  tha  index  of  refraction  a  »  1.33.  The  incident 
radiation  i3  along  the  symmetry  axis  (3/2  *  0)  and  4>.  *  3).  The  wave  number 
kQ  -  l  urn"1. 

In  these  thrae  figures,  the  second  order  produces  at  most  a  231  correction 
to  the  first  order  and  improves  agreement  with  the  exact  results.  The 
agreement  in  all  these  cases  is  excellent.  For  the  largest  spheroid  chosen, 
Figure  3.  the  quantity  3  *  k0Jl(m-L)  *  L  where  l  13  the  3pheroidal  parameter, 
t  *  /a^-b^.  Convergence  will  be  faster  for  S  <  L  and  it  will  be  slower  for  S 
greater  than  l.  A  rough  limit  on  the  useful  convergence  of  the  expansion  is 
3*2.  It  should  be  noted  that  this  criterion  for  convergence  applies  only  to 
the  Shifrin  technique  of  starting  from  the  electrostatic  limit. 


5.  COMPARISON  OF  THE  FINITE  CYLINDER  AND  SPHEROID 

The  amplitudes  Ij.  and  I^  ara  evaluated  using  the  electrostatic 
polarization  matrix  elements  of  the  short  cylinders.  These  results  for  the 
finite  dielectric  cyLinder  are  plotted  in  Figures  2  through  4  of  Reference  4  on 
which  the  values  for  the  spheroids  appear.  In  each  case,  the  aspect  catlos  and 
volumes  of  the  cylinder  and  3pheroid  are  the  same.  For  these  cases  the 
direction  of  incidence  is  along  the  symmetry  axis,  and  we  note  that  the 
cylinder  and  spheroid  give  remarkably  similar  scattering  patterns  for  tha 
aspect  ratios  1/2  and  l.  On  the  other  hand,  the  cylinder  with  aspect  ratio  5 
scatters  more  radiation  at  large  scattering  angles  than  Its  spheroid 
counterpart.  The  index  of  refraction  is  1.33  to  compare  with  tha  results  of 
Reference  3  (as  i3  the  selection  of  scattering  amplitudes  shown  in  these 
figures).  The  cylindrical  results  are  shown  as  dashed  lines  when  they  are 
different  from  the  spheroidal  results  (solid  lines).  The  open  circles  shown  in 
these  figures  are  exact  results.  The  wavelength  of  incident  radiation  is  2irym. 


Ia  Figure  2  of  Reference  4,  an  oblate  spheroid  and  cylinder  of  aspect 

ratio  1/2  are  shown,  each  with  a  volume  3.22  Urn  .  In  the  perturbative 

calculation  the  second-order  correction  is  generally  about  28?  for  both 

targets.  tfe  sea  that  the  scattering  pattern  for  the  amplitudes  I  and  I  are 
nearly  the  3ame.  1  2 

Figure  3  gf  Reference  4  compares  the  cylinder  with  a  sphere  each  with 

voLome  4.19ym  .  Here,  the  3acond-ocder  correction  is  generally  about  22?,  and 
the  results  for  the  cylinder  differ  noticeably  with  the  sphere  at  the  larger 
scattering  angles. 

Figures  4  and  5  of  Reference  4  show  the  respective  amplitudes  1^  and  l£ 
for  a  prolate  3pheroid  and  cylinder  of  aspect  ratio  5.  The  volume  of  each 
target  is  4.3lym-*.  The  second  order  correction  here  is  about  25?.  In  these 
cases,  the  cylinder  scatters  light  quite  differently  at  the  larger  scattering 
angles.  A  resonance  at  about  100°  appears  for  the  spheroid  in  the  amplitude  I 
but  does  not  appear  at  all  in  the  corresponding  cylinder.  2 

Figure  3  of  Reference  4  depicts  the  amplitude  for  incidence 

perpendicular  to  the  axis  of  symmetry  (g/2  »  90°,  <j>  *  9°).  Two  targets  are 

shown,  one  with  an  aspect  ratio  of  2  and  volume  l.52um-*,  and  the  other  has  an 
aspect  ratio  of  5  and  volume  4.3lya^.  The  second-order  correction  for  the 
first  is  about  20?  and  about  28?  for  the  longer  target.  We  note  that  the 
results  for  the  cylinder  differ  appreciably  from  the  corresponding  spheroid 
only  for  the  larger  aspect  ratio. 


5. 


ANGULAR  SCATTERING  DISTRIBUTIONS  BY  LONG  COPPER  AND  BRASS 
CYLINDERS:  EXPERIMENT  AND  THEORY 


Experimental  measurements  of  EM  radiation  scattered  by  long  copper  and 
brass  cylinders  ware  performed  in  the  IR  spectra  range  (X  *  l.Obym).  The 
cylinders  were  oriented  essentially  normal  to  the  scattering  plane.  Thi3  work 
was  done  at  Falrleigh  Dickinson  University  by  Tomasalli,  Moeller  and  Colosi. 
We  have  taken  the  results  of  the  measurements  and  compared  them  with  the  theory 
for  infinite  tilted  cylinders  modified  for  relatively  large  indices  of 
refraction.  The  far  field  scattering  of  infinite  tilted  cylinders  is  given  by 
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The  first  index  in  Iy  refers  to  the  polarization  of  the  incidence  light 
relative  to  the  incident  plane,  and  the  second-polarization  of  the  scattered 
light  raLative  to  the  scattering  plane.  For  detailed  definitions  of  the 
incident  and  the  scattering  planes,  3ae  Reference  9.  0*  is  the  scattering 
angle  and  bnj  can  be  reduced  to: 


AJ  (a)  +  [BJ'(a)  -  m2CJ  (a)][-  BH'(a)  -  CH  (a)] 
n _ n _  n _ n _ n _ 

AH  (a)  +  [BH'(a)  -  m2CH  (ct)][-  BH'(a)  -  CH  (a)] 
n  n  n  n  n 


(23) 


where 


nl 


a 

\  j 


h2(m2-l)2H  (a)J^(B) 
n  n 


B  =  *J  Jn(6) 


h  =  sin<|>  ;  l  =  cos<j>  ;  j  =  (mZ-sin2ij>)^2  ;  g  =  ^  a 


c  =  i2  jr  (g)  ; 


Jn(a)  and  Hn(a)  are  the  Bessel  and  Hankel  functions  of  order  n;  a  »  2-n-acos<f>/X  is 
the  "tilted"  size  parameter  (a  is  the  cylinder  radius,  X  the  incident 
wavelength  and  is  the  tilting  angle);  m  i3  the  complex  refractive  index.  The 
infinite  sums  in  equations  (20)-(22)  can  be  truncated  a  few  terms  after  the 
order  exceeds  the  "tilted"  size  parameter,  3ince  Jn(a)->0  for  a>p.  (n,  a 
Large).  For  example,  it  is  generally  agreed  that  for  a2  >  L  the  number  of 
terms  in  equations  (20)-(22)  is  of  the  order  of  N^l^ct+'S.  Therefore,  when  m  is 
large  (i.e.,  m  >  LO),  the  argument  ia  Jn(g)  satisfies  for  all  bnj  in  (20), 
g  »  n.  Or 


Jn(B)  * 


[cos(g  2)3 [1  +  0(|g|"1)] 


(24) 


For  copper  and  brass  cylinders  m  *  12-601  and  5.3-29i  respectively  (see  below). 
It  follows  that  for  any  size  parameter  a  *  2iracos<{)/X  larger  than  1,  Im(B)  be¬ 
comes  so  large  that  Jn(g)  cannot  be  calculated  since  it  contains  the  term 
exp  {lm(g)}.  Hence,  in  order  to  calculate  the  scattering  functions  for  those 
materials,  the  expression  for  bnI  has  to  be  modified: 

dividing  Equation  (23)  by  32  we  get 


J  («)  +  [Jl(c0  -  m2  £  J  (a) ] [ -  H’(a)  -  §  H  (a)] 
BZn  n  Bn  n  Bn 


nl  \  H  (a)  +  [H’(c0  -  m2  |  H  (a)  ]  [-  H'(a)  -  §  H  (a)] 
n  n  d  n  non 

o 


(25) 


The  ratio  G/ 3  is  the  only  term  containing  J'(g)  and  J^(g): 


Aj'(B) 
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For  large  values  of  g,  J^(8)  can  be  approximated  Co 

1/2 

^  -  -  (If)  •*>(»  -  T  -  3 

*=>  |  A#  ij  (Im(B)  <  0) 


We  note  that  for  m -►  «,  bni  reduces  to  the  known  ratio  bnj  **■  J  (a)/En(a). 
Similar  expressions  can  be  derived  for  AwXI  *  '^ea  the  cylinder  axis  is  not 
perpendicular  to  Che  measurement  plane  (see  Section  12),  the  scattering 
intensities  calculated  by  Equations  (29)-(22)  will  have  their  maximal 
intensities  within  the  envelope  of  a  cone  farmed  around  the  cylinder  axis.  In 
the  event  that  the  tilt  angle  between  such  a  cylinder  and  the  Incident 
direction  is  <j)  *  0,  the  scattering  plane  will  form  an  angle  $]_,  relative  to  the 
measurement  plane. 

The  comparison  between  theory  and  experiment  for  4>]_  «  0  then  becomes 
dependent  on  the  opening  angle  of  the  light  source  and  its  cross  sectional 
variation  of  Intensity  (usually  gaussian).  Here,  we  discuss  quantitatively  the 
ca3e  of  cf>i  »  0  (copper),  and  qualitatively  the  scattering  from  the  brass 
cylinder  for  which  <J) ^  /  0  but  £  i°. 


7.  OPTICAL  CONSTANTS  OF  TARGET  MATERIALS 

The  optical  constants  of  several  metals  from  infrared  to  far-infrared 
wavelengths  have  recently  been  tabulated  by  Ordal,  et  all®  Both  n  and  k  are 
rapidly  increasing  functions  of  wavelength  in  the  Infrared  region.  For  pure 
copper  at  10.6um,  representative  values  are  n  ■  12  and  k  ■  60.  The  brass  target 
wire  used  has  a  Cu/Zn  ratio  of  79/30.  For  metals  at  low  frequencies,  the 
optical  constants  can  be  approximated  by 


n  m  k  2  A/{ 2ue  o)  (2/) 

o 


where  p  is  the  static  resistivity,  co  is  the  angular  frequency,  and  £Q  i3  the 
permittivity  of  free  space.  Since  the  resistivities  of  metals  and  alloys  are 
easily  found,  we  used  the  above  approximation  together  with  the  tabulated  data 
for  p  to  obtain  the  optical  constants.  The  results  for  brass  at  lO.S^m  is 
n  *  5.3  and  k  ■  29. 

The  experimental  data  were  taken  for  wires  of  various  3ize3.  Figures  2a 
and  3a  of  Reference  6  represent  the  angular  scattering  results  from  copper  of 
diameter  242ym  and  brass  of  diameter  I50ym,  respectively. 


3.  RESULTS  AND  DISCUSSIONS 


The  theoretical  calculations  for  the  copper  wire  (a  »  1.2 1m)  show  that  the 
general  angular  scattering  as  well  as  the  accurate  angular  values  of  the  maxima 
and  minima  are  predictable  by  means  of  the  infinite  theory.  This  is  clearly 
shown  in  comparing  the  experimental  with  the  theoretical  results  in  Figure  2  of 
Reference  5.  The  accurate  match  shows  that  the  experiment  was  performed  in 
such  a  way  that  the  scattering  plane  and  the  measurement  plane  coincided  to 
permit  the  angular  measurements  in  the  3ame  relative  units.  It  is  important  to 
note  that  the  general  behavior  is  very  sensitive  to  the  size  parameter  and 
therefore  provides  an  accurate  method  for  the  determination  of  the  size  of  the 
scatterer.  However,  the  scattering  intensities  of  the  different  metallic 
cylinders  are  relatively  insensitive  to  changes  in  the  refractive  index.  This 
suggests  that  many  properties  can  be  approximately  derived  by  inserting  m  »  l 
in  the  scattering  equations. 

In  the  case  of  the  brass  cylinder  (the  experimental  results  are  given  in 
Figure  3a  and  the  theoretical  results  in  Figure  3b  of  Reference  5)  the 
scattering  angles  for  which  maxima  and  minima  occur  can  be  predicted,  but  the 
general  behavior  of  the  experiment  (a  decreasing  envelope)  i3  different  than 
the  theoretically  predicted  increasing  envelope.  However,  the  theoretical 
prediction  was  calculated  for  a  scattering  plane  that  is  perpendicular  to  the 
cylinder  axis.  In  the  event  that  the  wire  is  not  perpendicular  to  the 
measurement  plana,  having  an  orientation  angle  of  (ir/2  -  <J>),  the  scattering 
direction  will  be  tilted  relative  to  the  measurement  plane  with  a  varying  tilt 
angle  against  the  scattering  angle. 

As  can  be  seen  in  Figure  4  of  Reference  6,  this  varying  angle  reaches  its 
maximum  at  0  *  tt/2  where  it  equals  <p.  More  specifically,  denoting  the  measured 
scattering  angle  by  <f>  and  the  varying  tilt  angle  <j>'  ,  we  get 


tan  <J>'  »  3in  9  tan  4> 


(23) 


The  larger  the  angle  <j>'  is,  the  lass  is  the  overlap  between  the  scattering 
plana  field  at  view  (FOV)  and  the  measuring  plane  FOV  resulting  qualitatively 
in  a  general  decrease  of  the  scattering  intensity  with  intensity  increasing  in 
the  interval  0  <_  9  <_  ir/2. 

tfhen  <b  *  3,  the  scattering  angle  must  also  be  modified.  Denoting  the 
scattering  angle  in  a  plane  perpendicular  to  the  cylinder  by  <J>'  ,  the  following 
relations  nold  (see  Figure  4  of  Reference  5) 


sin  <j> '  *  sin  9  ’  sin  <p 


sin  9 ' 


cos  ' 
cos  $ 


sin  9 


As  stated  in  the  section  describing  the  experiment,  the  orientation  of  the 
wire  is  uncertain  to  within  3°.  Thus  < p  %  3°,  and  as  £  6,  the  deviation 
between  9'  and  9  is  less  than  the  accuracy  in  the  measurement  of  the  scattering 
angle  (»  4.3°).  Therefore,  thi3  effect  was  not  included  in  the  calculation. 
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9 .  MULTIPLE  SCATTERING 


As  a  first  step  to  a  mors  complete  description  of  the  multiple  scattering 
process  we  have  undertaken  to  describe  the  double  scattering  process  in  as 
complete  detail  as  possible.  The  geometry  of  the  double  scattering  process 
from  non-spherical  particles  is  particularly  complex,  however,  once  the 
geometry  is  worked  out,  the  transformations  necessary  to  consider  higher  order 
scattering  are  of  the  same  form  as  the  double  scattering  case. 

tfe  have  treated  in  detail  the  double  scattering  process  from  infinitely 
long  randomly  oriented  cylinders.  We  have  used  the  exact  infinite  cylinder 
solutions  so  that  the  cylinders  may  be  of  arbitrary  diameter  and  arbitrary 
index  of  refraction  (both  real  and  complex). 

The  geometry  for  calculating  the  backscatterlng  of  a  lidar  signal  from  an 
aerosol  layer  composed  of  randomly  oriented  long  cylinders  is  similar  to  the 
geometry  described  in  Reference  1.  Two  FOVs  are  of  interest  -  the  laser  FOV 
containing  the  illuminated  volume  and  hence  the  first  scatterer  and  the 
receiver's  FOV  containing  the  second  3catterer.  As  described  in  Reference  11, 
the  second  FOV  is  assumed  to  be  larger  or  equal  to  the  laser  FOV  (see  Figure 
1). 

The  aerosol  layer  i3  assumed  to  have  a  constant  number  density,  its  base 
being  1300m  aoove  the  lidar.  Since  the  laser  pulse  width  is  of  the  order  of 
100  nanoseconds  the  spatial  resolution  is  approximately  15m.  Consequently,  the 
optical  depth  3teps  were  based  on  the  product  of  the  extinction  cross-section 
of  a  typical  scatterer  averaged  over  all  possible  orientations,  and  the  number 
density  multiplied  by  15m.  Reference  ll  deals  with  a  given  size  distribution 
of  spherical  scatterers.  Here,  we  deal  with  long  cylinders  randomly  oriented 
in  space  but  having  one  given  size. 

In  order  to  calculate  the  sum  of  all  possible  double  scattering  events 
resulting  in  a  simultaneous  signal  in  the  receiver,  a  statistical  approach  i3 
used. 


Two  randomly  chosen  scattering  positions  inside  the  cloud  [called  A  (the 
lower  point)  and  3  (the  higher  point)  in  the  following  analysis]  are  determined 
by  the  computer.  The  first  has  the  constraint  of  being  within  the  telescope 
FOV.  An  additional  constraint  requires  that  the  sum  of  the  distances  from  the 
source  to  A,  A  to  B,  and  3  to  the  receiver  be  constant  to  within  the  resolution 
of  the  system  for  simultaneous  reception.  In  each  position,  it  is  assumed  that 
all  possible  orientations  can  be  found;  however,  each  orientation  in  the  first 


location  corresponds  to  one  scattering  angle  pointing  towards  the  second 
location.  This  is  due  to  the  fact  that  the  light  scattered  by  a  unit  length  of 
the  cylinder  is  limited  to  the  surface  of  a  cone  based  on  the  assumption 
mentioned  above  that  the  scatterers  are  effectively  infinite  cylinders  for 
aspect  ratios  greater  than  about  100  (sea  Reference  3  and  the  discussion  of 
Reference  5).  The  scattering  events  A  and  3  are  shown  in  Figure  2  with  these 
two  points  chosen  at  random  (with  B  higher  than  A  by  definition).  The  incident 
direction  is  ,  the  first-scattered  direction  is  ,  and  the  final  scattered 
direction  is  .  The  reference  plane  is  formed  by  z  and  z^g  .  The  incident 
wave  denoted  ,  and  its  direction  of  polarization  makes  an  angle,  a ,  (relative 
to  the  reference  plane)  as  seen  in  Figure  2.  The  angle  between  the  zQ  axis  and 
AB  i3  2$0 ,  and  we  note  from  figure  2  that 


m 


*  jN’  * 
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Figure  1 

Two  OF  THE  FAMILY  OF  CONES  THAT  CORRESPOND  TO  A 
GIVEN  PAIR  OF  SCATTERERS , 


with  xQ  ini  RAg  shown  La  Figure  2,  with  0  <_  4>0  VJe  will  next  show  that 
the  tilt  angles  and  <j>g  of  the  two  cylinders  are  the  only  orientational 
degrees  of  freedom.  Note  that  only  those  cones  whose  tilt  angles  i.  satisfy 
4) A  >  4>Q  2 -if1  contribute  to  the  ioubLe  scattering,  as  shown  In  Figure  3. 


10.  0RISNTATI0N  CONSTRAINTS 


We  define  the  tilt  angle  of  a  cylinder  to  be  the  angle  between  the  axis  of 
the  cylinder  and  the  incident  direction  of  the  light  (the  C^-axis)  -  these 
angles  are  denoted  4>  and  for  the  two  cylinders  as  3hown  In  Figure  2,  and 
0  ^  (j>  <_  7T.  A  B 

First,  we  note  that  far  the  section  of  a  cone  shown  In  Figure  2  with  apex 
angle  4>a»  24>q  the  angle  between  the  extreme  elements  of  the  section,  and  0A  the 
polar  angle  at  the  base  of  the  section 


sin  4>a  sin(0A/2) 


sin  4> 


Xhi3  relation  holds  for  the  first  conic  section  at  the  point  A  in  FLgure  2. 
The  corresponding  relatLon  for  the  conic  section  at  the  point  B  Is 


sin  4>n  s^n  (80/2)  ■  cos  4>  (31) 

D  D  O 

In  the  following  discussloa,  we  use  a  set  of  orientation  angles  a,  yA*  snd  Yg- 
All  of  these  are  measured  by  a  counterclockwise  rotation  from  the  reference 
pLane.  In  Figure  4a,  the  unit  vector  A  follows  the  direction  of  the  ^axis  of 
the  cylinder  at  point  A  with  its  origin  at  paint  A.  The  unit  vector  B  follows 
the  direction  of  axis  Using  the  frame  of  reference  xD,  yQ,  zQ,  shown  In 

the  fLgure,  we  have  the  following  representations: 


A  =  sin  <p.  cos  y,  i  +  sin  4>.  sin  y.  j  +  cos  d>.  k 

A  A  o  A  A  o  A  o 


B  *  sin  24>  i  +  cos  2<p  k 
o  o  o  o 


The  scalar  product  then  yields 


cot  4>»  ■  sin  24)  cos  y.  +  cot  4>.  cos  24> 


(32) 


cot  4.  =  cot  $  COS  Y. 

A  o  A 


>  tt/2 . 


(Note  that  0  <_  y  <_  tt  ,  and  since  cot  $0  >  0,  when  y  >_  t/2  than  iJ:A  >_  tt/2. 
Therefore,  of  the  three  orientation  angles  for  the  first  scattering  event  (i 
the  tiLt  angLe  far  the  cylinder,  0  the  scattering  angle,  and  yA  the  azimuthal 
angle  of  the  cylinder  axis  as  shown  in  Figure  4),  there  is  only  one  that  is 
independent  by  the  relationships  of  Equations  (30)  and  (32).  Similar  results 
hold  for  the  second  cylinder. 

Two  additional  angles  are  required  to  describe  the  double  scattering  - 
angle  between  the  first  scattering  plane  and  the  Incidence  plane  for  the  second 
scattering,  and  the  angle  between  the  second  scattering  plane  and  the  direction 
of  initial  polarization.  These  angles  are  easily  determined  in  terms  of  the 
orientation  angles  already  discussed  by  first  determining  the  orientations  of 
these  scattering  planes  relative  to  the  reference  plane. 

In  Figure  4b,  the  angle  ya  i3  the  angle  between  the  xD  axis  in  the 
reference  plane  and  the  first  incidence  plane  (the  first  incidence  plane  i3 
formed  by  the  incident  direction  and  the  cylinder  axis).  We  have  already 
represented  the  unit  vector  A  that  i3  along  the  axi3  of  the  first  cylinder  in 


terms  of  the  frame  xr 


angle  yA  can  be  determined  from  the  ratio 


zQ,  and  it  follows  from  thi3  that 


azimuthal 


cot  Y. 


A  /A 
x  y 
o  1  o 


The  angle  <5A  i3  the  corresponding  azimuthal  angle  in  the  frame  xAB,  yQ, 
zAB,  and  hence  satisfies  a  similar  expression  in  this  new  frame  ( 3ea  Figure 
4b).  These  two  frames  of  reference  are  related  by  a  counterclockwise  rotation 
of  xQ ,  yQ ,  zQ  through  the  angle  2<pQ  about  the  yD  axis: 
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It  follows  that 


cot  6  a  =  axab/\ 


cos  2o  cot  y.  “  cot  sin  2j>  esc  y. 
o  A  A  o  A 
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A  3Lnil3r  result  holds  for  the  orientation  of  the  second  scattering  plane 
(formed  by  the  axi3  of  the  second  cylinder  and  the  final  scattering  direction), 
and  this  situation  is  shown  in  the  Figures  5a  and  5b. 

The  next  section  discusses  the  scattering  intensities.  We  will  need  the 
following  angles  to  complete  that  discussion. 


Angle  from  the  incident  polarization  direction  to  the  incident  plane 
rotated  about  z  in  a  counterclockwise  direction  *  yA  -  a 
Angle  between  the  scattering  plane  to  the  second  incidence  plane  rotated 
about  zAB  =»  YB  -  <SA  a  Yn  -  Ya, 

Angle  between  the  second  scattering  plane  and  the  incident 
polarization  direction  rotated  about  -  zQ  ■  -  a  -  <Sg  =  -  a  -  Yg. 

A  study  of  Figures  4  and  5  should  clarify  these  relationships. 


U.  00U3LE  SCATTERING  INTENSITY 

The  transmitted  beam  and  the  double  scattering  intensities  are  expressed 
below  in  terms  of  the  modified  Stokes  vector.  The  transmitted  beam  is  a 
linearly  polarized  laser  beam  with  intensity 


I  -  (I  0,  0,  0) . 
o  o 


The  intensity  of  the  lidar  in  the  plana  of  polarization  i3  denoted  IQ,  and  the 
double  scattering  intensity  components  Xs 1 1  ,  Is^>  are  expressed  relative  to  I0. 
In  ocder  to  calculate  the  double  scattering  intensities  in  the  geometry 
discussed  above,  the  cloud  volumes  contained  in  the  F3V  angles  of  the  laser  and 
the  receiver  are  divided  into  n  identical  3ub-volumes.  The  two  scattering 
events  occur  in  a  given  cloud  layer  of  depth  L,  situated  at  a  height  H  (above 
the  lidar  system).  The  double  scattering  intensity  is  calculated  by  adding  the 
intensities  scattered  by  all  possible  pairs  of  the  two  sub- volumes  (one  in  the 
transmitter  and  one  in  the  collector  FOV),  for  which  the  total  path  length  of 
the  light  equals  the  average  path  to  the  top  of  the  cloud  layec  and  back  to  the 
lidar  system  (3ea  Figure  1). 

The  power  scattered  at  an  angle  9^  emerging  from  a  sub-volume  AVA,  is 
described  by  the  vector  I^(9A»  <J>A)  related  to  the  incident  light  vector  IQ  by 
the  scattering  matrix  for  cylinders  P  (9,6),  and  a  rotational  matrix 
50vo  - 1>  as  follows:  ^ 


ava 

IA(0A’  V  “  Pcy(9A’  V  <X^A  *  a)Io  ZT  Gcy-'AB  '  "0A' 

rab 


(35) 


with  a  Che  extinction  coefficient.  The  scattering  matrix  P  (9,6)  is 
cy  5  cv  ’  * 


Figure  5 


General  geometry  for  double  scattering  of  lidar  radiation 

FROM  A  DENSE  CLOUD  OF  PARTICLES. 
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rtia  matrix  elements  ace  related  to  the  cylinder  scattering  amplitudes  Alt  \2, 
A^  which  are  functions  of  the  scattering  and  tilt  angles  (0A,  <j>A) : 
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Ine  pathlength  for  the  Lidar  system 

to  the 

first  scattering  sub- 

denoted  RQA,  and  the  pathlength  between  the  scattering  sub-volumes  is  denoted 


AB’ 


Then,  using  the  discussion  of  the  last  section  on  geometry,  we  can 


express  the  intensity  of  the  scattered  Light  after  the  double  scattering  events 


(36) 


F  AV  AV 
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where  aCy  is  the  extinction  coefficient  of  the  randomly  oriented  cylinders,  F  is 
the  lidar  constant,  is  the  distance  between  the  two  scattering  sub-volumes 

^Ai»  aa(*  aBoi  ace  tlia  ai3Caace  between  the  second  scattering  volume  AVg^  and 
the  lidar  system. 

Note  that  >  AR/10,  with  AR  the  spatial  lidar  resolution,  therefore  no 

singular  points  are  present.  The  condition  R^g^  <  AR/IG  is  regarded  as  single 
scattering.  The  pathlength  within  the  cloud  that  does  not  include  R^m  13 
denoted  Rcf  The  double  scattering  geometry  is  discussed  in  the  previous 
section.  In  aquation  (36)  the  summation  i3  only  over  the  pairs  of  sub-volumes 
obeying 


a«wli  *■  aA6>  f  a5<Jl  *  l2s 


(37) 


Here,  RQA  is  the  pathlength  from  the  lidar  system  to  the  Etrst  scattering 
volume,  and  a  is  the  height  of  the  top  of  the  cloud  layer  above  the  cloud  base. 
Since  the  size  of  the  sub-volumes  should  be  as  small  as  possible  in  order  to 
reduce  averaging  effects,  the  summation  over  the  number  K  of  all  pair  (A,  B) 
satisfying  equation  (3/)  requires  extensive  computer  time  for  a  case  of  double 
scattering.  For  this  reason,  a  large  sample  of  g(«  K)  pairs  of  points  are 
chosen  la  six  random  steps  (which  determine  the  coordinates  of  the  scattering 
points  A  and  B) .  Each  selection  places  the  point  A  in  the  first  sub-volume  and 
3  in  the  second,  and  each  pair  of  points  replaces  a  pair  of  sub-volumes.  As 
the  total  number  of  points  N  (not  necessarily  satisfying  equation  (37))  is 
luiown,  the  value  of  X  can  be  calculated  by  means  of  the  random  program.  Each 
pair  of  random  points  is  checked  to  satisfy  equation  (37)  in  the  following  way: 
the  ratio  G  *  a/g,  where  a  is  the  number  of  pairs  contributing  to  the  double 
scattering  process,  is  recorded  and  therefore  K  ■  NG. 

In  addition  to  the  3ix  random  steps  which  determine  the  position  of 
scattering  cylinders,  yA  and  yg  are  randomly  chosen  from  which  the  scattering 
and  polarization  angles  are  calculated  as  discussed  in  the  last  section.  In 
order  to  get  the  double  scattering  Intensity  profile  from  the  whole  cloud, 
calculations  were  performed  for  several  depths  in  steps  corresponding  to  the 
spatial  resolution  of  a  conventional  lidar  system.  In  each  step,  the  whole 
random  process  described  above  was  evaluated  changing  the  value  of  Rg  in 
equation  (37). 

Since  extinction  and  multiple  scattering  Intensities  are  strongly 
dependent  on  the  number  density  of  the  scatters,  calculations  were  performed 
for  various  relative  densities,  and  the  results  are  presented  as  intensity 
profiles  as  a  function  of  the  penetration  depth  R  in  the  cloud  for  various 
number  densities  (Figures  6  and  7). 
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Figure  7 

Double  scattering  intensities  as  a  function  of  penetration  depth 
FOR  NUMBER  DENSITIES  FROM  1100  TO  1500  PARTICLES  /  CM?  . 
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L2.  DISCUSSION 


The  double  scattering  effect  of  randomly  oriented  long  cylinders  was 
calculated  for  number  densities  varying  from  5x10^  to  10^  and  3hown  la  Figure 
5  where  0.2  £  T  £  0.4,  sad  for  number  densities  10^  to  1.5x10^  as  shown  In 
Figure  7,  where  4  <  t  <  0.5.  In  general,  the  behavior  of  the  double  scattering 
profile  as  a  function  of  the  penetration  depth  resembles  the  scattering  by 
dense  clouds  of  spherical  particles.  .4.3  can  be  seen  in  Figure  S,  the  double 
scattering  effect,  for  relatively  low  number  densities  is  to  produce  a  maximum 
scattering  layer  which  approaches  the  cloud  base  with  an  increased  value  of 
number  density.  Single  and  single  plane  double  scattering  intensities  are 
shown  in  Figure  3  for  the  number  densities  500  to  1000  particle/cm^.  The 
double  scattering  readouts  are  obtained  from  the  results  shown  in  Figure  5. 

It  is  important  to  emphasize  that  the  geometrical  approach  discussed  here 
is  the  first  such  analysis  to  appear  in  literature.  Moreover,  it  can  be 
applied  to  other  non-spherical  particles  for  which  the  scattering  matrices  are 
expressed  la  terms  of  a  symmetry  axis.  For  example,  this  treatment  is  easily 
extended  to  finite  cylinders  by  using  the  matrices  developed  in  Reference  3. 
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_  Single  scattering  intensities 

-  Single  plus  double  scattering  intensities 
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